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We continue our exploration of whether the flyby anomalies can be explained by scatter- 
ing of spacecraft nucleons from dark matter gravitationally bound to the earth, with the 
addition of data from five new flybys to that from the original six. We continue to use our 
model in which inelastic and elastic scatterers populate shells generated by the precession of 
circular orbits with normals tilted with respect to the earth's axis. With 11 data points and 
8 parameters in the model, a statistically meaningful fit is obtained with a chi-squared of 
2.7. We give plots of the anomalous acceleration along the spacecraft trajectory, and the cu- 
, mulative velocity change, for the five flybys which exhibit a significant nonzero anomaly. We 

• ■ also discuss implications of the fit for dark matter-nucleon cross sections, give the prediction 

of our fit for the anomaly to be expected from the future Juno fiyby, and give predictions of 
our fit for fiyby orbit orientation changes. In addition we give formulas for estimating the 
flyby temperature increase caused by dark matter inelastic scattering, and for the fraction of 
flyby nucleons undergoing such scatters. Finally, for circular satellite orbits, we give a table 
of predicted secular changes in orbit radius. These are much too large to be reasonable - 
comparing with data for COBE and GP-B supplied to us by Edward Wright (after the first 
version of this paper was posted), we find that our model predicts changes in orbit radius 
that arc too large by many orders of magnitude. So the model studied here is ruled out. We 
conclude that further modeling of the flyby anomalies must simultaneously attempt to flt 
constraints coming from satellite orbits. 



I. INTRODUCTION 

In this paper we follow up our earlier investigations [1], [2] of the anomalous geocentric frame 
orbital energy changes that are observed during earth flybys of various spacecraft, as reported 
by Anderson et al. [5]. Some flybys show energy decreases, and others energy increases, with 
the largest anomalous velocity changes of order 1 part in 10^. While the possibility that these 
anomalies are artifacts of the orbital fitting method used in [5] is still being actively explored, 
there is also a chance that they may represent new physics. 
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In [1] we explored the possibility that the flyby anomalies result from scattering of spacecraft 
nucleons from dark matter particles in orbit around the earth, with the observed velocity decreases 
arising from elastic scattering, and the observed velocity increases arising from exothermic inelastic 
scattering, which can impart an energy impulse to a spacecraft nucleon. Many constraints on 
this hypothesis were analyzed in [1], with the conclusion that the dark matter scenario is not 
currently ruled out, but requires dark matter to be non-self-annihilating, with the dark matter 
scattering cross section on nucleons much larger, and the dark matter mass much lighter, than 
usually assumed. These constraints on the dark matter scattering cross section and mass will be 
revisited, and somewhat weakened, in the present paper. 

In [2] we constructed a model for the spatial and velocity distribution functions for dark matter 
particles in earth orbit, based on assuming two populations of dark matter particles, one of which 
scatters on nucleons elastically, and the other of which scatters inelastically, with each population 
taken to fill a shell-like distribution of orbits generated by the precession of a tilted circular orbit 
around the earth's rotation axis. Fit of this model to the data were presented, but since the 
model has 8 parameters, and there were only 6 data points in the original flyby set, the statistical 
significance of the fits was doubtful. After [2] was written, Anderson and Campbell [6] furnished 
us with data for 5 additional fiybys, all of which show no velocity anomaly. In the present paper, 
we add these 5 fiybys to the original 6 and fit the model of [2] to all of the data points. We find 
that with 11 data points and 8 parameters, we get a statistically meaningful fit with a chi-squared 
of 2.7. In this fit, the dark matter shells lie much closer to earth than they did in the fits exhibited 
in [2]. 

The formulas defining the model were given in Sec. II of [2] (with further detail in the Ap- 
pendices of the arXiv version), and will not be repeated here. So reading the earlier paper is a 
necessary prerequisite for fully understanding the model analyzed in this paper. Our focus here 
is on presenting the orbital parameters of the new fiybys and the results of the fit to the full 11 
flyby set. In addition to giving the fits to the total velocity anomaly, we give the contributions of 
the shells of elastic and inelastic scatterers to this total. We also plot the anomalous acceleration 
along the spacecraft trajectory and the cumulative velocity change, for the five fiybys which show 
significant nonzero velocity anomalies, give constraints on dark matter-nucleon cross sections im- 
plied by our fit, and give formulas for estimating the fiyby temperature change and the fraction of 
nucleons that undergo exothermic dark matter scatters. For the future Juno flyby, our 11 flyby fit 
predicts an anomaly of 11.6 mm/s. 

We conclude by giving predictions of our model for fiyby orbit orientation anomalies, and for 
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secular changes in orbit radius of satellites in circular orbits. The latter, obtained from a program 
that can treat non-circular orbits as well, shows that the predicted orbit radius changes for COBE 
and GP-B are orders of magnitude larger than observational data supplied to us by Edward Wright 
after the first version of this paper was posted to arXiv: see the Added note to Section VII. Hence 
the model studied here is ruled out, and any further fitting of more general forms of the dark 
matter model, or of other types of models, to the flyby anomalies, must simultaneously attempt 
to fit constraints coming from closed satellite orbits. This greatly complicates the analysis, and 
we shall pursue this further only if the future Juno flyby, which has orbital parameters similar to 
GLL-I, confirms the presence of a kinetic energy anomaly. 



II. FLYBY ORBITAL PLANE PARAMETERS 



As in [2], it will be convenient for to carry out all flyby orbit calculations in the flyby orbital 
plane. Let Xq, Uo, Zq be a Cartesian axis system, with Zg normal to the flyby orbital plane. The 
flyby orbit can then be written in parametric form as 

Xo (t) =r {t) cos 9o{t) , yg{t) = r (t) sin 9o (t) , 



1 + ecos 6'o(t)' 1 + e ' 

1 + e l + ecost^o(t) 



(1) 



The scale parameter p, the eccentricity e, the velocity at closest approach to earth V/, the radius 
at closest approach Rf, and the velocity at infinity Fqo are given in Table I for the six original 
fiybys, and in Table II for the five new ones and the future Juno flyby, together with the polar 
angle / and azimuthal angle a of the earth's north pole with respect to the Xo,yo,Zo coordinate 
system. Formulas for obtaining Rf,p,e, and a from the data presented in [5] and [6] are given in 
[2]. 
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TABLE I: Orbital parameters for the original six flybys 





GLL-I 


GLL-II 


NEAR 


Cassini 


Rosetta 


Messenger 


Vf (km/s) 


13.740 


14.080 


12.739 


19.026 


10.517 


10.389 


Rf (km) 


7,334 


6,674 


6,911 


7,544 


8,332 


8,715 


Voo (km/s) 


8.949 


8.877 


6.851 


16.010 


3.863 


4.056 


e 


2.474 


2.320 


1.814 


5.851 


1.312 


1.360 


p (km) 


25,480 


22,160 


19,450 


51,690 


19,260 


20,570 


/ (deg) 


142.9 


138.7 


108.0 


25.4 


144.9 


133.1 


a (deg) 


-45.1 


-147.4 


-55.1 


-158.4 


-53.1 


0.0 



III. MODEL PARAMETERS AND FIT TO 11 FLYBYS 

The model of [2] postulates two shells of dark matter gravitationally bound to earth. An inelastic 
scatterer shell is the locus obtained by starting with a circular orbit with normal tilted at angle 
with respect to the earth's rotation axis, and smearing this orbit into a uniform shell by rotation 
around the earth's rotation axis to represent the effect of quadrupole moment-induced precession. 
The radial density in the shell is assumed to have a Gaussian profile proportional to e~^^~^^^^/^i , 
and the effective shell density times scattering cross section on fiyby nucleons is denoted by pi. 
This gives four parameters for the inelastic shell, the tilt angle i/^i, the central shell radius Ri, 
the shell width Dj, and the effective density pi. For the elastic scatterer shell (which could arise 
as daughter products from interactions of the inelastic scatterers with nucleons) we assume an 
analogous structure, introducing four more parameters, again a tilt angle ipe, a central radius Rg, 
a shell width D^, and an effective density p^. A more detailed description of these parameters is 
given in Section 2.6 of [2] (Section IIF of the arXiv version). 



TABLE IL Orbital parameters for the five new flybys, and the future Juno flyby 





Rosetta II 


Rosetta III 


EPOXI I 


EPOXI II 


EPOXI III 


Juno 


Vf (km/s) 


12.49 


13.34 


6.94 


5.29 


5.73 


15.024 


Rf (km) 


11,656.6 


8,860.62 


21,985.02 


49,786.17 


36,775.84 


6871.48 


Voo (km/s) 


9.36 


9.38 


3.45 


3.46 


3.34 


10.474 


e 


3.562 


2.956 


1.656 


2.495 


2.029 


2.891 


p (km) 


53,177.89 


35,051.16 


58,402.95 


174,016.7 


111,402.94 


26738.31 


/ (deg) 


115.08 


155.63 


19.83 


92.81 


103.07 


48.23 


a (deg) 


175.35 


110.04 


158.28 


139.34 


-142.29 


140.01 
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We note that the closest approach to earth of the flybys in our fit (flyby GLL-II) is i?/ = 6, 674 
km, corresponding to an altitude H = 303 km, and thus the fits would be unaltered if the shell 
radial profile were changed to a cut-off Gaussian e~^^~^*)^/^< 0(r — 6,674km), with 9{r) the usual 
step function. Hence the model cannot be used to make predictions for orbits that lie below an 
altitude of 303 km, which includes the majority of low earth orbit satellites. Predictions of the 
model for earth orbiting satellites in circular orbits lying above the cut-off will be given below. 

As in [2] , fitting this model to the data was carried out by minimizing a least squares likelihood 
function x^, defined as 

11 

= ^{6vk;th - 5vk-Af/al.^ , (2) 
fc=i 

where k indexes the six original flybys discussed in [5] and the five new flybys reported in [6], 
where the 5vk-t\i are the theoretical values of the velocity anomalies computed from our model, 
the 5vfe;A are the observed values for these anomalies reported in [5] and [6], and the a^.^x are 
the corresponding estimated errors in these anomalies. Since the quoted a^.^x values contain both 
systematic and statistical components, a least squares likelihood function is not a true statistical chi 
square function, but having a quadratic form is very convenient for the following reason. Because 
the theoretical values (5ffe;th are linear in the dark matter density times cross section parameters 

Pi,e-i 

Svk-th = PiSvk-i + PeSvk,e , (3) 

with 6vk-i^e tbe respective contributions from the inelastic and elastic scatterers computed with 
Pi^e = 1) the likelihood function is a positive semi-definite quadratic form in these two parameters. 
Hence for fixed values of the other six parameters ipi^e, Ri,e-i -f^i,e5 the minimization of with 
respect to the parameters pi^^ can be accomplished algebraically by solving a pair of linear equations 
in the two variables pi^e-, as described in Section 3 of [2] (Section HI of the arXiv version). 

As a consequence, it is only necessary to search numerically a six parameter space. We performed 
the fit by doing a survey of the six parameter space on a coarse mesh to find a good starting point for 
the minimization program Minuit, then doing a finer mesh survey centered on the minimum found 
by Minuit, giving an improved starting point for another application of Minuit, and so forth, until 
this process converged. The results are given in Tables HI through VH. Tables HI and IV give the 
input experimental values of the velocity anomalies and their estimated errors, together with the 
fits obtained by our search program (using an adaptive integration method). The chi-squared value 
at the minimum was 2.7, and the values of the eight parameters at the chi-squared minimum are 
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TABLE III: Experimental values and 11 fiyby fit 





GLL-I 


GLL-II 


NEAR 


Cassini 


Rosctta 


Messenger 


Sy-A (mm/s) 


4.48 


-4.60 


13.46 


-1.02 


1.80 


0.02 




0.3 


1.0 


0.01 


1.0 


0.03 


0.01 


6v;th (mm/s) 


4.09 


-4.66 


13.5 


-0.807 


1.80 


0.01 



TABLE IV: Experimental values and 11 flyby fit continued, and prediction for the future Juno flyby 





Rosetta II 


Rosetta III 


EPOXI I 


EPOXI II 


EPOXI III 


Juno 




(mm/s) 


0.0 


0.0 


0.0 


0.0 


0.0 






o-A 


0.1 


0.1 


0.1 


0.1 


0.1 




Sv;th 


(mm/s) 


0.0 


0.006 


0.0 


0.0 


0.0 


11.6 



given in Table V. Using these parameters, we also recalculated the theoretical velocity anomalies 
using a 200,001 point trapezoidal integration, as well as printing out the separate contributions of 
the elastic and inelastic shells to the total theoretical value, as given in Tables VI and VIL We see 
that the total theoretical anomaly is approximately independent of the integration method, and 
that the elastic and inelastic shells both make important contributions in our model. 

In about two years time, the Juno flyby should give an additional data point; the prediction of 
our fit for this flyby is a large anomaly of 11.6 mm/s, as also shown in Table IV. 



TABLE V: Parameter values for 11 flyby fit 



10'^ X Pi (km) 


10^ X pe (km) 


ipi (rad) 


V'e (rad) 


Ri (km) 


Di (km) 


Re (km) 


De (km) 


0.398 


0.272 


2.79 


0.0603 


7561 


2038 


12526 


1668 
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TABLE VI: 11 flyby fit (tot) and contributions from inelastic (i) and elastic (c) shells 





GLL-I 


GLL-II 


NEAR 


Cassini 


Rosctta 


Messenger 


Sv;th tot (nim/s) 


3.97 


-4.65 


13.4 


-0.755 


1.76 


0.010 


Sv;th i (mm/s) 


8.11 


7.63 


13.5 


14.3 


8.98 


0.010 


6y-th e (mm/s) 


-4.14 


-12.3 


-0.03 


-15.1 


-7.22 


0.00 



IV. PLOTS OF THE INELASTIC AND ELASTIC SHELL ACCELERATION ALONG 
THE TRAJECTORY, AND OF THE CUMULATIVE VELOCITY ANOMALY 

In this section we describe sample plots, given at the end of the paper, of the spacecraft down- 
track acceleration from the inelastic and elastic dark matter shells (i.e., the inelastic and elastic 
shell accelerations projected along the spacecraft trajectory) and the cumulative asymptotic veloc- 
ity anomaly 6V(x, (obtained by integrating the total downtrack acceleration with respect to time, 
dividing by the velocity at infinity, and multiplying by 10^ to get an answer in mm/s rather than 
in km/s).^ Plots are given for the five fiybys GLL-I, GLL-II, NEAR, Cassini, and Rosetta which 
had significant nonzero velocity anomalies. For NEAR, we also plot the downtrack, crosstrack, and 
normal total accelerations, defined respectively as the inner product of the total acceleration with 

1 /2 

(dxo/dt, dyo/dt, 0)/v, {—dyo/dt, dxo/dt, 0)/v, and (0,0,1), with v = (^{dxo/dt)'^ + (dyo/dt)'^^ 
the downtrack velocity. We also give plots of the cumulative asymptotic velocity anomaly and the 
downtrack acceleration from the inelastic and elastic shells for the future Juno flyby, for which a 
large anomaly is expected based on our fit. Because of the integrable singularities in the accelera- 
tions arising from the Jacobian that is discussed in detail in [2] , we have plotted the accelerations on 
a semi-log rather than a linear plot; on a linear plot, only a spike at the Jacobian peaks appears in 
any detail. These plots may facilitate comparison of the dark matter model with other anomalous 



TABLE VII: 11 flyby fit (tot) and contributions from inelastic (i) and elastic (e) shells - continued 





Rosetta II 


Rosctta III 


EPOXI I 


EPOXI II 


EPOXI III 


Sv;th tot (mm/s) 


0.0 


0.008 


0.0 


0.0 


0.0 


Sv;th i (mm/s) 


0.0 


5.14 


0.0 


0.0 


0.0 


Sv;th e (mm/s) 


0.0 


-5.13 


0.0 


0.0 


0.0 



^ Note that the cumulative asymptotic velocity anomalies exhibited are inferred from the corresponding kinetic 
energy anomahes in the asymptotic region where the potential energy is negligible; these graphs do not give 
the changes in the along-orbit velocity in sub-asymptotic regions where the potential energy must be taken into 
account. See Appendix A for a detailed discussion. 
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force models devised to explain the flyby anomalies; on request by email to adler@ias.edu, we can 
also furnish the numerical tables from which these plots were made. 

V. PREDICTIONS FOR FLYBY ORIENTATION ANOMALIES 

Using the formulas given in Appendix A, we can calculate the anomalies in the angular mo- 
mentum vector per unit mass L, and the Laplace-Runge-Lenz vector per unit mass squared A, 
produced by passage of a flyby through the dark matter shells. Together with the change in the 
energy, or equivalently the change in the along-trajcctory asymptotic velocity, on which we have 
based the fit to our model, these additional anomalies completely characterize the final asymptotic 
state of the flyby. The results of this calculation are given in Table VIII. 



TABLE VIII: Anomalies in 5L{oo) in km^/s, and in 5A{oo) in km^/s^, with components given on the x,y, 
and normal z axes of the flyby plane. 



flyby 




SLy 


6L, 




SAy 


SA, 


GLL-I 


0.13E-01 


-0.23E-01 


0.12E-01 


0.47E+00 


-0.17E+00 


-0.12E+00 


GLL-II 


-0.32E-01 


-0.27E-01 


-0.47E-01 


-0.77E+00 


0.65E+00 


0.31E+00 


NEAR 


0.19E-01 


-0.31E-01 


0.52E-01 


0.13E+01 


0.20E+00 


-0.16E+00 


Cassini 


-0.29E-01 


-0.22E-01 


0.26E-01 


0.29E+00 


-0.42E+00 


0.47E+00 


Rosctta 


0.96E-02 


-0.25E-01 


0.20E-02 


0.13E+00 


-0.16E+00 


-0.57E-01 


Messenger 


O.OOE+00 


O.OOE+00 


O.OOE+00 


O.OOE+00 


O.OOE+00 


O.OOE+00 


Rosetta II 


O.OOE+00 


O.OOE+00 


O.OOE+00 


O.OOE+00 


O.OOE+00 


O.OOE+00 


Rosetta III 


-0.14E-02 


0.14E-01 


-0.31E-02 


-0.25E-01 


-0.12E+00 


0.14E-01 


EPOXI I 


-0.40E-26 


0.56E-24 


0.12E-22 


0.16E-21 


-0.16E-23 


0.17E-25 


EPOXI II 


O.OOE+00 


O.OOE+00 


O.OOE+00 


O.OOE+00 


O.OOE+00 


O.OOE+00 


EPOXI III 


O.OOE+00 


O.OOE+00 


O.OOE+00 


O.OOE+00 


O.OOE+00 


O.OOE+00 


Juno 


-0.14E-01 


0.19E-01 


0.64E-01 


0.18E+01 


0.36E+00 


0.16E+00 
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VI. FORMULAS FOR THE FLYBY TEMPERATURE CHANGE, AND THE 
FRACTION OF SCATTERED NUCLEONS 

In addition to flyby velocity changes arising from the average over the scattering cross section 
of the coUision-induced velocity change, there will be spacecraft temperature increases arising from 
the mean squared fluctuation of the collision-induced velocity change. These are estimated in 
Appendix B (which simplifies an earlier account given in [4]), giving for inelastic scattering the 
formula 



_ _ 0.13°K Jl^mm'c^ 

^ ' I (cos 0)1 MeV ' ^ ^ 

with m!^ the dark matter secondary mass, with Am = — m'2 the dark matter exothermic energy 
release, and with 6 the nucleon scattering angle. Additionally, taking the ratio of the flyby velocity 
change ~ 1 cm/s ~ 0.3 x 10~^°c to Eq. (B14) of Appendix B for the nucleon velocity change in a 
single inelastic scatter, we see that the fraction F of flyby nucleons undergoing inelastic scatters is 
of order 

„ 0.2x10-10 mi 

F ~ . (5) 

|(cos6')| y/Amm^ ' ^ ' 

with mi the nucleon mass. These formulas can be used to place constraints on the model, when 
upper bounds on flyby temperature changes, and on possible radiation damage to flyby electronic 
components, are available. 



VII. PREDICTIONS FOR ORBIT RADIUS CHANGES OF EARTH-ORBITING 

SATELLITES 

The formulas of Eq. (1) also describe closed elliptical orbits when e < 1, with the semi-major 
axis a given by a = p/(l — e^). So the same program used to evaluate the flyby anomalies can also 
be used to evaluate the secular change 5a /a for satellites in closed orbits, induced by scattering on 
the dark matter shells. To get the secular change over a single orbit, we take the limits of the 9o 
integration for the energy change 5E as — tt < 9o <t^, and use the formula 

5a 2a5E 

Sample results of this calculation, for circular orbits (e=0) and with the result re-scaled to give the 
secular change over a year, are given in Table IX. We see that the secular changes depend strongly 
on the orbit, and can range from very small to large. For polar orbits, they are relatively small. 
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and for orbits beyond a radius of 20,000 km they are negligible. This table shows that it will be 
important to impose closed orbit constraints on the model. The data needed for this (for satellites 
beyond an altitude of 303 km) are the orbit semi-major axis a, the ellipticity e, and the earth axis 
polar angle / and azimuthal angle a on the satellite orbit plane, together with the observed value 
of, or bound on, the annual anomalous increment baja. 

Added note: The above paragraph is what appeared in the first version of this paper. As a 
result of our posting this to arXiv, Edward Wright sent us several emails with detailed data for 
COBE, GP-B and other satellites. The secular changes in orbit radius for these satellites are orders 
of magnitude smaller than predicted by Table IX. For example, for COBE, which is in a near-polar 
orbit with an inclination of 99 degrees, and a=7278 km, the observed dln{a)/dt = — 0.00009/yr, 
compared with a prediction of order 0.02/yr from Table IX. And for GP-B, which is in an almost 
exactly polar orbit with inclination of 90.007 degrees, and a=7027.4 km, quoting from Wright's 
email, "from 2004.7 to 2004.8 GP-B was in a "drag-free" mode, where the spacecraft tracked the 
motion of one (of) the rotors that was shielded from atmospheric drag and radiation pressure. I have 
attached a plot of the orbital rate of GP-B vs date. The drag-free epoch is quite obvious.... during 
the drag-free epoch \dln{a)/dt\ is clearly less than 0.000005/yr, while your model predicts about 
0.028/yr." We agree with Wright that this data rules out the model analyzed in this paper. Any 
further modeling of the fiyby anomalies must, as an essential component, include closed orbit 
satellite constraints in the fitting procedure. Our perspective is that this should wait until after 
the Juno flyby, which has an orbit similar enough to that of GLL-I (compare Figs. 1-3 for GLL-I 
with Figs. 19-21 for Juno) to have the potential, independent of modeling, to confirm or rule out 
the presence of a kinetic energy anomaly. If the Juno flyby can be tracked through earth approach, 
and shows an anomaly, analysis of the tracking data using the methods of Appendix A would give 
important information about the radius range where the anomaly arises. 



In [2] we gave an estimate for the total mass Me and Mj in the elastic and inelastic dark matter 



VIII. DISCUSSION 



shells. 



Me c^ATT^/'^peDemi/aei 




(7) 
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where in Mj we have allowed for the possibility that Am is much smaller than m'2 (rather than 
Am ~ as assumed in [2]). Here a^i is the threshold cross section for the scattering of the elastic 
dark matter population on nucleons, B-^^^i is the coefficient of the cos 9 term in the near threshold 
cross section for exothermic scattering of the inelastic dark matter population on nucleons (see 
Eq. (7) of [2]), and mi, m^^ and Am are as defined in Section VI and Appendix B. Using these 
estimates, and the upper bound [3] of 4 x lO^^M® ~ 1.4 x 10^'^GcV/c^ on the mass of dark matter 
in orbit around the earth between the 12,300 km radius of the LAGEOS satellite orbit and the 
moon's orbit, one can get lower bounds on a^i and Anei- For example, referring to Table V, from 
PeDe = 4.537 km^ we get 

(Tel > 2.2 X 10-^^cm^ , (8) 
while from piDi = 0.0000811 km^, taking into account the fact that only a fraction 

1 r°° 

/ dzexp{-z^) ~ 0.50 X 10"^ (9) 

(with L = (12, 300 - 7561)/2038 = 2.325) of the inelastic dark matter sheh lies above the LAGEOS 
orbit, we correspondingly get 

Anei-^Am/m^ > 1.4 x lO'^^cm^ . (10) 

These bounds arc consistent with the cross section range arrived at from various constraints in 
[1]. In comparing these bounds with cross section limits derived from dark matter direct detection 
experiments, two points should be kept in mind. The first is that if the dark matter mass is much 
lower than a GeV, the recoils from elastic scattering on nucleons of galactic halo dark matter will be 
below the threshold for detection in direct searches, and so no useful constraints on the elastic dark 
matter nucleon scattering cross section are obtained. The second is that with two species of dark 
matter, as in our model, the fraction /j of the galactic halo consisting of inelastically scattering dark 
matter can be much smaller than unity. In this case a direct detection cross section bound (Jmax, 
obtained assuming a single species of dark matter in the galactic halo, implies a significantly larger 
bound <TMAx//i for the dark matter-nucleon scattering cross section of the inelastic component. 

As opposed to the fits given in [2], the current fit to 11 flybys places the elastic and inelastic dark 
matter shells much closer to earth, well within the orbits of geostationary and Global Positioning 
System satellites. So the effect on these satellites should be minimal, but it now becomes important 
to examine the implications of our model for satellites that are in low or medium earth orbit, as 
initiated in Section VII and Table IX. (See in this regard the Added note in Section VII.) It will 
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also be important to evaluate implications of a low earth orbit exothermic scattering dark matter 
shell for earth's heat balance. We pose these as significant questions to be addressed in future 
investigations. 
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Appendix A: Calculation of energy, angular momentum, and Laplace— Runge—Lenz vector 
anomalies, and their relation to local velocity and position anomalies 

When a perturbation of force per unit mass 6F is applied to a flyby spacecraft, the total change 
in energy per unit mass at time t is given by the work integral 

6E{t)= f du^^= f duvo{u)-5Fiu)= f du^^^-6Fiu) , (Al) 

which when divided by the velocity at infinity and multiplied by 10^ is the quantity plotted as 
the "cumulative asymptotic velocity anomaly SV^o" in the figures. Taking a first variation of the 
energy per unit mass (where we write ro{t) = \xo{t)\ and /x = GM®), 

we get 

ft 



6E(t)= [ du^^^-6F(u) 
J-oo du 

=SVo[t) ■ Vo{t) + --3 



(A3) 
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This formula gives the relation between the cumulative energy anomaly at time t, and the position 
and velocity anomalies of the orbit Sxo{t) and 6vo{t) at time t. At time t = oo, the potential energy 
contribution vanishes, giving a formula relating the work integral taken from t = — ootoi = oo 
to the asymptotic along-track velocity anomaly SV^o, which as this derivation makes clear is an 
asymptotic flyby energy anomaly. 

In addition to the energy, a Kepler orbit has conserved angular momentum vector per unit mass 

L, 

L = Xo{t) X Vo{t) , (A4) 
and conserved Laplace-Runge-Lenz vector per unit mass squared A, 

A = Vo{t)xL-^^ . (A5) 
Taking the time derivatives of these quantities, substituting the perturbed equation of motion 

and then integrating with respect to time, and equating these to first variations of the conserved 
quantities, we get formulas 

SL{t) = / duXo{u) X 5F{u) 

J —oo 

=5Xo{t) X Vo{t) + Xo{t) X 5Vo{t) , 

(A7) 



and 



5A{t) = I du^5F{u) ■ Vo{u)xo{u) — 6F{u) ■ Xo{u)vo{u) — Xo{u) ■ Vo{'(J')SF{u)] 

J —CXI 

=2Vo{t) ■ 5Vo{t)Xo{t) + Vo{tf5Xo{t) - SVoit) ■ Xo{t)Vo{t) - Vo{t) ■ 5Xo{t)Voit) 
^ ^ ll5Xo{t) fJ,Xo{t)Xo{t) ■ 5Xo{t) 

-v^{t) . x^{t)Sv^{t) - + . 

(A8) 

Setting t = oo in these formulas gives, in terms of integrals over the unperturbed orbit of the force 
perturbation SF, an expression for the asymptotic orbit orientation perturbations 5L{oo) and 
5^(00). These in turn can be calculated using the above formulas from the observed asymptotic 
perturbations 6xo{oc) and 5vo{oo) obtained from tracking the flyby position and velocity, 

SL{oo) = 5xo{oo) X 'yo(oo) + Xo(too) x Svo{oo) , (A9) 
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and 

SA(oo) =2vo(oo) ■ 5voioo)xo{oo) + Vo{oo)'^6xo{oo) — dvo(oo) ■ Xo{oo)vo{oo) 
—Vo{oo) ■ Sxo{oo)vo{oo) — Voioo) ■ Xo{oo)Svo{oo) 

(AlO) 

As a point of consistency, we note that the same values of 6L{oo) and dA{oo) are obtained irrespec- 
tive of where on the asymptotic outgoing trajectory these formulas are evaluated. If after a first 
evaluation, they are then evaluated at a time T later on, the quantities Xo(oo) and 6x0(00) are aug- 
mented respectively by Vo{oo)T and 5vo{oo)T. One can then check that the additions proportional 
to T cancel out of the equations for 6L{oo) and 6A{oo). 

Appendix B: Temperature change eirising from velocity fluctuations 

In [1] we considered the velocity change when a spacecraft nucleon of mass mi 2± IGeV and 
initial velocity ui scatters from a dark matter particle of mass m2 and initial velocity U2, into an 
outgoing nucleon of mass mi and velocity vi, and an outgoing secondary dark matter particle of 
mass m'2 = m2 — Am and velocity V2 ■ (In the elastic scattering case, one has m2 = 1712 and 
Am = 0.) Under the assumption that both initial particles are nonrelativistic, so that |ui| << 
c,\u2\ « c, a straightforward calculation shows that the outgoing nucleon velocity is given by 

VI = ■ h WVout • (Bl) 

mi + 1712 

Here it; > is given^ by taking the square root of 

2 m2m2 ,^ ^ .2 , Am m'2 [ 2 {miui+m2U2f 1 ^^^s 

(mi -|- m2j(mi + m2j mi (mi -|- m2j L (mi -|- m2j(mi + m2j-l 

and 'Oout is a kinematically free unit vector. Denoting by 9 the angle between -Oout and the entrance 
channel center of mass nucleon velocity ui — (miui +m2U2)/{mi +m2) = m2{ui — U2)/{mi +m2), 
and assuming that the center of mass scattering amplitude is a function f{9) only of this polar 
angle, the average over scattering angles of the outgoing nucleon velocity is given by 

miUi + m2U2 , , ^.ui-U2 , 

{vi} = ■ \-w{cose}-^ ^ , (B3) 

mi -|- m2 1^1 — U2\ 



^ The notation t was used in [1] for what we here term w; the change in notation avoids confusion with use of t for 
time. 
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with (cos 9) given by 

de sine coseif (0)1'^ , , 

Subtracting u\ from Eq. (B3) gives the formula for the average velocity change used in [1] and [2] 
to calculate the flyby velocity change, 

\ov\) = = V w Leo's, B)-—, rri ■ (B5) 

mi + m2 1^1 — U2I 

However, in addition to contributing to an average change in the outgoing nucleon velocity, 
dark matter scattering will give rise to fluctuations in this velocity, which have a mean square 
magnitude given by 

, 2 

This fluctuating velocity leads to an average temperature increase of the nucleon, per single scat- 
tering, of order 



((^1 - (^1))') = ^'(kut - if^osB) ^} = v?{\ - {cosef) . (B6) 



with kB the Boltzmann constant. In analogy with the treatment of the velocity change 5vi in [1], 
to calculate dT/dt, the time rate of change of temperature of the spacecraft resulting from dark 
matter scatters, one multiplies the temperature change in a single scatter (ST) by the number 
of scatters per unit time. This latter is given by the flux \ui — U2\, times the scattering cross 
section cr, times the dark matter spatial and velocity distribution p{x,U2\ Integrating out the 
dark matter velocity, one thus gets for dT/dt at the point x{t) on the spacecraft trajectory with 
velocity ui = dx{t)/dt, 

dT/dt = J d^U2{ST)\ui-U2\ap{x,U2) ■ (B8) 

Integrating from ti to tf we get for the temperature change resulting from dark matter collisions 
over the corresponding interval of the spacecraft trajectory , 



Tf — Ti= dt d'^U2{ST)\ui — U2\ap[x,U2) ■ 
Jti J 



(B9) 



In the elastic scattering case, with Am = 0, 1712 = 1712, the formula of Eq. (B2) simplifies to 
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In the inelastic CclSGj clS long as Am/m2 >> \ui — U2\ /c , Eq. (B2) is well approximated by 

2 / 2Am 1712 \ 2 

W ~ 7 ^ c . (Bll) 

\ mi(mi + I 

Since ui and U2 arc typically of order 10 km s^^, the temperature change in the inelastic case, per 
unit scattering cross section times angular factors, is larger than that in the elastic case by a factor 
~ c^/|nip 10^, and the restriction on Am is the weak condition Am/m2 » 10~^. 

We proceed now to estimate Eq. (B9) by using Eqs. (B7), (BIO), and (Bll), and making the 
approximations that the dark matter mass m2 is much smaller than the nucleon mass mi , and that 
(cos 6') 2 in Eq. (B7) can be neglected relative to 1. In the elastic case, Eq. (4) of [1] tells us that 
the magnitude of the velocity change in a single collision is of order 

|(M)| - — |ui -■U2I . (B12) 
mi 

Taking the ratio of the single collision temperature change to the single collision velocity change, 
and multiplying by the flyby total velocity change ~ 10~^|mi|, we get as an estimate of the total 
temperature change 



in agreement with [1]. 

In the inelastic case, we must take into account the kinematic structure of an exothermic inelastic 
differential cross section. In the inelastic case, Eq. (5) of [1] tells us that the magnitude of the 
velocity change in a single collision is of order 



\{6vi)\^^^^^^c\{cos9)\ . (B14) 
mi 



Writing the inelastic differential cross section near threshold in the form 

dcr ^inel , _i „ 3 



we have 



A;-^ + Bineii-cos6' + (B15) 



(C0S6') ~5inel/(Anelfe~^) , 



(B16) 



with k the entrance channel momentum 



K= \Ul — U2\ — m2\Ul — U2\ ■ (-Dl'j 

77ll + 7712 
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Again taking the ratio of the single colhsion temperature change to the single collision velocity 
change, and multiplying by the flyby total velocity change ~ lO^^lt/il, with l^il ~ lOkm/s, we get 
as an estimate of the total temperature change in the inelastic case 

^ ^ ST _ fi,^ , 10-6 1 \ui\ r To 0.13°K JAmmLc^ 

^ \{5vi)\ ' ' 3A;b |(cos6')| c V 2 |(cos6l)| MeV ^ ' 

So wc see that the exothermic inelastic scattering temperature rise is substantially bigger than 
that from elastic scattering, as already anticipated in the remarks following Eq. (Bll) above. In 
particular, for the temperature change to be bounded, say, by 10°K, Eq. (B18) implies that 

^J/^mm'^c^ < 77|(cos6»)|MeV , (B19) 



which is compatible with y'Amm^ in the McV range. Wc note however, that if Am/m2 were much 
smaller than unity, then the dark matter mass m2 could be much larger than an MeV. 
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TABLE IX: Change 5a/ a for a circular orbit in a year, versus orbit radius and earth orbit polar angle. The 
first column is the orbit radius in km, with the remaining columns giving the annual 5a /a as a function 
of the earth polar axis angle / listed at the top of the table. For / = 0, the orbits circulate in the same 
direction as the earth's rotation, and for I = 180 the orbits circulate in the opposite sense as the earth's 
rotation. For circular orbits, the results are independent of the earth axis azimuthal angle a. 



radius 


0.0 


22.5 


45.0 


67.5 


90.0 


112.5 


135.0 


157.5 


180.0 


6878. 


0.16E+00 


0.21E+00 


0.73E-01 


0.44E-01 


0.29E-01 


0.20E-01 


0.12E-01 


0.70E-02 


0.51E-02 


7378. 


0.15E+00 


0.20E+00 


0.71E-01 


0.42E-01 


0.28E-01 


0.19E-01 


0.12E-01 


0.67E-02 


0.47E-02 


7878. 


0.13E+00 


0.18E+00 


0.61E-01 


0.37E-01 


0.24E-01 


0.16E-01 


O.lOE-01 


0.56E-02 


0.31E-02 


8378. 


O.lOE+00 


0.14E+00 


0.47E-01 


0.28E-01 


0.19E-01 


0.12E-01 


0.77E-02 


0.38E-02 


-0.13E-02 


8878. 


0.70E-01 


0.93E-01 


0.32E-01 


0.19E-01 


0.13E-01 


0.81E-02 


0.44E-02 


0.67E-03 


-0.14E-01 


9378. 


0.43E-01 


0.57E-01 


0.20E-01 


O.llE-01 


0.69E-02 


0.36E-02 


0.20E-03 


-0.52E-02 


-0.46E-01 


9878. 


0.24E-01 


0.31E-01 


O.lOE-01 


0.52E-02 


0.18E-02 


-0.15E-02 


-0.61E-02 


-0.17E-01 


-0.12E+00 


10378. 


0.12E-01 


0.15E-01 


0.41E-02 


0.42E-03 


-0.32E-02 


-0.79E-02 


-0.16E-01 


-0.37E-01 


-0.25E+00 


10878. 


0.50E-02 


0.61E-02 


0.23E-03 


-0.35E-02 


-0.85E-02 


-0.16E-01 


-0.29E-01 


-0.65E-01 


-0.44E+00 


11378. 


0.19E-02 


0.18E-02 


-0.22E-02 


-0.67E-02 


-0.14E-01 


-0.24E-01 


-0.43E-01 


-0.97E-01 


-0.65E+00 


11878. 


0.65E-03 


-0.71E-04 


-0.35E-02 


-0.88E-02 


-0.17E-01 


-0.30E-01 


-0.54E-01 


-0.12E+00 


-0.80E+00 


12378. 


0.18E-03 


-0.73E-03 


-0.39E-02 


-0.93E-02 


-0.18E-01 


-0.31E-01 


-0.56E-01 


-0.12E+00 


-0.84E+00 


12878. 


0.36E-04 


-0.80E-03 


-0.35E-02 


-0.81E-02 


-0.16E-01 


-0.27E-01 


-G.49E-01 


-O.llE+00 


-0.73E+00 


13378. 


-O.lOE-05 


-0.62E-03 


-0.25E-02 


-0.60E-02 


-O.llE-01 


-0.20E-01 


-0.36E-01 


-0.80E-01 


-0.53E+00 


13878. 


-0.61E-05 


-0.39E-03 


-0.16E-02 


-0.37E-02 


-0.70E-02 


-0.12E-01 


-G.22E-01 


-0.49E-01 


-0.33E+00 


14378. 


-0.41E-05 


-0.20E-03 


-0.80E-03 


-0.19E-02 


-0.36E-02 


-0.63E-02 


-O.llE-01 


-0.25E-01 


-0.17E+00 


14878. 


-0.19E-05 


-0.87E-04 


-0.35E-03 


-0.81E-03 


-0.16E-02 


-0.27E-02 


-0.49E-02 


-O.llE-01 


-0.73E-01 


15378. 


-0.71E-06 


-0.31E-04 


-0.13E-03 


-0.29E-03 


-0.56E-03 


-0.98E-03 


-0.18E-02 


-0.39E-02 


-0.26E-01 


15878. 


-0.22E-06 


-0.95E-05 


-0.38E-04 


-0.89E-04 


-0.17E-03 


-0.30E-03 


-0.53E-03 


-0.12E-02 


-0.80E-02 


16378. 


-0.55E-07 


-0.24E-05 


-G.96E-05 


-0.22E-04 


-0.43E-04 


-0.75E-04 


-0.13E-03 


-0.30E-03 


-0.20E-02 


16878. 


-0.12E-07 


-0.51E-06 


-0.20E-05 


-0.48E-05 


-0.91E-05 


-0.16E-04 


-0.29E-04 


-0.64E-04 


-0.43E-03 


17378. 


-0.21E-08 


-0.91E-07 


-0.36E-06 


-0.85E-06 


-0.16E-05 


-0.28E-05 


-0.51E-05 


-O.llE-04 


-0.76E-04 


17878. 


-0.31E-09 


-0.14E-07 


-0.54E-07 


-0.13E-06 


-0.24E-06 


-0.42E-06 


-0.75E-06 


-0.17E-05 


-O.llE-04 


18378. 


-0.39E-10 


-0.17E-08 


-0.67E-08 


-0.16E-07 


-0.30E-07 


-0.52E-07 


-0.94E-07 


-0.21E-06 


-0.14E-05 


18878. 


-0.40E-11 


-0.18E-09 


-0.70E-09 


-0.16E-08 


-0.31E-08 


-0.55E-08 


-0.98E-08 


-0.22E-07 


-0.15E-06 


19378. 


-0.35E-12 


-0.15E-10 


-0.61E-10 


-0.14E-09 


-0.27E-09 


-0.48E-09 


-0.85E-09 


-0.19E-08 


-0.13E-07 


19878. 


-0.26E-13 


-O.llE-11 


-0.44E-11 


-O.lOE-10 


-0.20E-10 


-0.35E-10 


-0.62E-10 


-0.14E-09 


-0.93E-09 


20378. 


-0.16E-14 


-0.68E-13 


-0.27E-12 


-0.63E-12 


-0.12E-11 


-0.21E-11 


-0.38E-11 


-0.85E-11 


-0.57E-10 
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acceleration from inelastic shell km/s^ 
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FIG. 2: GLL-I downtrack acceleration from inelastic shell km/s^ 
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cumulative velocity anomaly mm/s 
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FIG. 4: GLL-II cumulative asymptotic velocity anomaly SVoo in mm/s 
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FIG. 6: GLL-II downtrack acceleration from elastic shell km/s^ 
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cumulative velocity anomaly mm/s 




FIG. 7: NEAR cumulative asymptotic velocity anomaly SVoo in mm/s 




FIG. 8: NEAR downtrack acceleration from inelastic shell km/s^ 




FIG. 9: NEAR downtrack acceleration from elastic shell km/s^ 



total crosstrack acceleration km/s^ 




FIG. 11: NEAR total crosstrack acceleration km/s^ 



31 



cumulative velocity anomaly mm/s 




FIG. 13: Cassini cumulative asymptotic velocity anomaly 6Voo in mm/s 
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acceleration from elastic shell km/s^ 
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FIG. 15: Cassini downtrack acceleration from elastic shell km/s^ 
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FIG. 16: Rosetta cumulative asymptotic velocity anomaly dVoo in mm/s 



acceleration from inelastic shell km/s 




FIG. 17: Rosetta downtrack acceleration from inelastic shell km/s^ 
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FIG. 18: Rosetta downtrack acceleration from elastic shell km/s^ 
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cumulative velocity anomaly mm/s 




FIG. 19: Juno cumulative asymptotic velocity anomaly 6Voo in mm/s 




FIG. 21: Juno downtrack acceleration from elastic shell km/s^ 



